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Abstract 



Oh| The measurement of hadronic Higgs Boson branching ratios H— )• 

^ ' bb, H— )• gg for a hght Standard Model-hke Higgs boson produced at 

250 GeV centre of mass energy at the International Linear Colhder 
(ILC) is presented. The tools and techniques used for the analysis are 
^ ■ briefly discussed. 



1 Introduction 



m 

O I The measurement of Higgs boson branching ratios is one of the main features 

of the International Linear Collider (ILC) program [Ij . For Higgs masses be- 
low 140 GeV, hadronic branching ratios can be precisely measured at the 
ILC. The final states have significant rates and the micro-vertex detector al- 
lows for good flavour identification. This measurement of relative couplings 
<xi I of the Higgs boson to fermions will allow to confirm the prediction of the 

Higgs mechanism that they are proportional to fermion masses. The branch- 
ing fraction to b and c quarks is larger than the branching fraction to light 
quarks due to the mass. The Higgs decay to gluons in the Standard Model 
is mediated by heavy quark loops as shown in Figure [H The branching ratio 
to gluons is indirectly related to UR Yukawa coupling [2] and would probe 
the existence of new strongly interacting particles that couple to the Higgs 
and are too heavy to be produced directly. 

In this study we consider Higgs decays to bottom quarks and gluons in 
two- and four-jet configurations. For a mn = 120 GeV SM Higgs boson, we 
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Figure 1: Higgs decay to gluons mediated by a heavy quark loop. 

expect BR(H ^ bb) = 67.92% and BR(H ^ gg) = 7.06% [3]. These decay 
modes exercise the tagging of heavy and hght quarks [H |5]. 

2 Analysis 

The analysis strategy closely follows the one implemented in 

2.1 Event Generation and Detector Simulation 

In this study, the signal sample includes a Higgs boson produced in the 
Higgs-strahlung process, e"'"e~— )■ ZH. The mass of the Higgs is assumed to be 
120 GeV. Standard Model events (mainly WW, ZZ and qq pairs) and Higgs 
decays to other fermions other than gluons/bottom quarks are considered 
as backgrounds. Signal and background events are produced at the centre 
of mass energy y/s = 250 GeV, total integrated luminosity of 250 fb~^ and 
=F80% electron and ±30% positron polarization |4j. The choice of energy in 
this analysis maximizes the cross-section value for Higgs-strahlung. All 0, 
2 and 4 fermion final states were generated with the Whizard Monte Carlo 
Event Generator [B]. PYTHIA [7] was used for the final state QED and QCD 
parton showering, fragmentation and decay to provide final-state observable 
particles. Photons from beamstrahlung and initial state radiation were also 
included in the simulations. 
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Geant4 toolkit [HI E] was used to simulate detector response to generated 
events and SLIC [lU] provided access to the Monte Carlo events, the detector 
geometry and the output of the detector hits. 

2.2 Event Selection 

The event selection is identical to that done in [1] where further analysis 
details can be found. There are two channels studied in this analysis: the 2-jet 
neutrino channel (when the Z boson decays to neutrinos and the Higgs decays 
to gg/bb) and the 4-jet hadronic channel (Z decays to quark pairs and Higgs 
to gg/bb). The selection of events defining each channel is based on the visible 
energy and the number of leptons in the event. Hadronic jet reconstruction, 
in which events are forced into two or four jet configurations, is achieved 
by the DURHAM algorithm [11]. Identification of primary, secondary and 
tertiary vertices is performed by the topological vertex finder ZVTOP which 
is part of the vertexing package developed by the LCFI collaboration [T2] . 

After the channel classification, a neural network analysis is performed 
in order to discriminate the signal and background events. The surviving 
events are then used for the branching ratio calculation. 

2.2.1 Neutrino Channel 

The signal in this channel consists of two jets from the Higgs recoiling against 
two neutrinos from the Z boson decay. The missing mass is expected to be 
consistent with the Z mass and the invariant mass of the 2 jets consistent 
with the Higgs mass. The main backgrounds in this channel include WW 
pairs where 

• one W decays hadronically. 

• and the other W decays into a neutrino and lepton which escapes un- 
detected along the beampipe, 

ZZ pairs in which one Z decays hadronically and the other into neutrinos and 
qq pairs. 

For this channel, no leptons are accepted and the visible energy is required 
to be between 90 and 160 GeV. Leptons are defined as electrons or muons 
with minimum 15 GeV momentum. The following pre-selection cuts are 
applied to reduce the background: 

1. 20 < Pt (transverse momentum) of jet < 90 GeV. Most SM background 
events are softer than signal events. 
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2. ntracks > 4. Moie than 4 charged tracks for leptonic event rejection. 

3. — \og{ymin) < 0.8. Durham algorithm parameter which determines 
number of jets in events. It is used to reject fully hadronic WW and 
ZZ events. 

4. thrust < 0.95. Background events are more boosted and less spherical 
than signal events. 

5. cos(9thrust)< 0.98. Signal events occur more centrally in the detector 
than background events. 

6. 100° < angle between jets < 170°. 

7. 100 GeV < di-jet invariant mass < 140 GeV. The Higgs mass is ex- 
pected to be 120 GeV. 

8. Highest reconstructed photon energy < 10 GeV. Required to reject 
2-fermion events with large ISR. 

Tables [T] and [2] show the number of events before and after pre-selection 
cuts for bb and gg in the neutrino channel. 



Cuts 


SM background 


Higgs background 


Signal 


(i) Before Classification 


9275594683 


6048 


12187 


(0) After Classification 


45936973 


5248 


11580 


(1) 


18374789 


5053 


11243 


(2) 


17123140 


4255 


10609 


(3) 


6849256 


3976 


10196 


(4) 


685329 


3782 


8427 


(5) 


627113 


3562 


8027 


(6) 


576422 


3403 


7907 


(7) 


203292 


2786 


6801 


(8) 


109057 


2737 


6707 



Table 1: Number of bb events before channel classification, after channel 
classification and after pre-selection cuts in the neutrino mode. 



Events that survive the pre-selection cuts undergo a neural network (NN) 
analysis to help discriminate the signal and background. The input variables 
to the neural network include all pre-selection variables and also the jet 
flavour tagging outputs produced using the LCFI package. There are three 
types of jet flavour tag outputs: b-tag, c-tag and c-tag with b background 
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Cuts 


OA /T "U 1 J 

SiVl background 


TT* "U 1 J 

Higgs background 


O * 1 

Signal 


(i) Beiore ulassmcation 


9z75594o83 


17547 


992 


Alter Uiassitication 


459oo97o 


15820 


986 


(1) 


18374789 


15317 


979 


(2) 


17123140 


13896 


968 


(3) 


6849256 


13252 


920 


(4) 


685329 


11343 


865 


(5) 


627113 


10766 


823 


(6) 


576422 


10515 


795 


(7) 


203292 


8818 


769 


(8) 


109057 


8684 


759 



Table 2: Number of gg events before channel classification, after channel 
classification and after pre-selection cuts in the neutrino mode. 



only. Figure [2] and Figure [3] show distributions of the three LCFI flavour 
tags, 'b-tag', 'c-tag' and 'c-tag with b background only' for the leading b and 
gluon jets. For the bb/gg scenarios the signal is defined as only H— t- bb/gg 
events and the Higgs background is all other Higgs decays other than H— t- 
bb/gg. All histograms are normalized to 250 fb~^. 

The tagging of jets works well. B-jets have relatively long lifetimes and 
decay topology with up to two secondary vertices and this makes them very 
distinct. Gluonic jets are tagged as light quark jets. 

2.2.2 Hadronic Channel 

In the hadronic channel both the Higgs and Z bosons decay into two partons 
(quarks or gluons) resulting in a 4-jet configuration. For the signal, the 
invariant mass of two jets is required to be consistent with the Higgs mass 
and the other two jets should have a mass corresponding to the Z boson. To 
reduce combinatorial effects, kinematic fitting [H [13] is performed to identify 
jets from the Higgs boson and jets from the Z boson. Jet pairing is performed 
before kinematic fitting. For the 4-jet events we have 6 possible pairings of 
the jets and 3 possible associations of the 4 jets to the Z and H bosons. 
For the 6 possible pairings we calculate the invariant mass of each pair and 
compare with the masses of bosons. For each event we calculate 

d = {rriij - nizf + {mti - ninf (1) 

The pairing that minimizes d is chosen. The signal four-jet configuration 
has, on average, the leading and third jets coming from the Higgs boson and 
the second and last jets coming from the Z boson. 
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(a) 



(b) 



I total tagged bcs jet 1 | 




0.2 0.4 0.6 0.8 



(c) 

Figure 2: Distributions of flavour tagging variables for bb leading jet in the 
neutrino channel: (a) b-tag; (b) c-tag; (c) c-tag with b background only. Solid 
curves are SM background, dashed curves are Higgs background sample and 
filled histograms are the signal. 
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I total tagged bs jetT 



I total tagged cs jetT 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



(a) (b) 

total tagged bcs jet 1 | 

10' 
10' 
10» 
10' 



(c) 

Figure 3: Distributions of flavour tagging variables for gg leading jet in neu- 
trino channel: (a) b-tag; (b) c-tag; (c) c-tag with b background only. Solid 
curves are SM background, dashed curves are Higgs background sample and 
filled histograms are the signal. 
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The kinematic constraints used for the fitter are the jet four-momenta, 
centre of mass energy (250 GeV) and the invariant mass difference of two jet 
pairs. 

Classification of the hadronic channel requires that no leptons are present 
in the event and a minimum of 170 GeV of visible energy. Apart from the 
variables used in the neutrino channel, the other variables/cuts applied in 
the hadronic channel are the angle between Z boson jets and the invariant 
mass of jets coming from the Z boson. Table [3] shows the selection cuts in 
the hadronic channel. 



Cuts 






selection 


value 


1. 






number of charged tracks per jet 


> 4 


2. 






- \og{ymin) 


< 2.7 


3. 






thrust 


< 0.95 


4. 






COs{6thrust) 


< 0.96 


5. 


105° 


< 


angle between jet 1 and 3 


< 165° 


6. 


70° 


< 


angle between jet 2 and 4 


< 160° 


7. 


110 GeV 


< 


invariant mass of Higgs candidate after fit 


< 140 GeV 


8. 


80 GeV 


< 


invariant mass of Z candidate after fit 


< 110 GeV 


9. 






Highest reconstructed photon energy 


< 10 GeV 



Table 3: Selections for the four-jet analysis. 



Tables |4] and [5] show the number of events before and after pre-selection 
cuts for bb and gg in the neutrino channel. 



Cuts 


SM background 


Higgs background 


Signal 


(i) Before Classification 


9275594683 


17816 


35629 


(0) After Classification 


39398366 


11875 


30985 


(1) 


18601753 


8821 


24398 


(2) 


13921271 


7582 


20736 


(3) 


8737017 


6472 


18391 


(4) 


7943851 


5977 


17130 


(5) 


5871237 


5681 


16339 


(6) 


4898312 


5678 


16326 


(7) 


1917231 


5633 


16108 


(8) 


1561432 


5622 


16108 


(9) 


967312 


5486 


15805 



Table 4: Number of bb events before channel classification, after channel 
classification and after pre-selection cuts in the hadronic mode. 
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Cuts 


OA /T "U 1 J 

SiVl background 


TT" "U 1 J 

Higgs background 


* 1 

Signal 


(i) Beiore ulassmcation 


9z75594o83 


50453 


2992 


(Uj Alter uiassmcation 


39398366 


39894 


2965 


(1) 


106U1753 


30717 


2503 


(2) 


13921271 


26119 


2199 


(3) 


8737017 


22977 


1887 


(4) 


7943851 


21396 


1738 


(5) 


5871237 


20373 


1648 


(6) 


4898312 


20357 


1648 


(7) 


1917231 


20103 


1638 


(8) 


1561432 


20092 


1638 


(9) 


967312 


19860 


1611 



Table 5: Number of gg events before channel classification, after channel 
classification and after pre-selection cuts in the hadronic mode. 



The main backgrounds at this stage are fully hadronic WW and ZZ pairs, 
and 2-fermion pairs. As in the neutrino channel, events that survive the pre- 
selection cuts undergo a neural network analysis to help further discriminate 
the signal and background. The input variables to the neural network include 
all pre-selection variables and also the jet flavour tagging outputs produced 
using the LCFI package. Figures HI [5], [6] and [7] show distributions of the 
three LCFI flavour tags, 'b-tag', 'c-tag' and 'c-tag with b background only' 
for the highest (leading) and second highest energy b- and gluon jets. 

2.2.3 Di-jet Mass Resolution 

Jet energy measurements are a crucial part of Higgs studies and the mea- 
surement of the Higgs mass distribution width is important. Figures [8] and 
[9] show the reconstructed visible mass distribution of Higgs decays to bb, cc 
and gg in the neutrino and hadronic channels respectively. The histograms 
are fit with a single gaussian which is not always adequate but nevertheless 
allows qualitative analysis. In the neutrino channel the hh and cc systems 
have broader mass widths than the gluons because of extra neutrinos coming 
from b- abd c-hadron semi-leptonic decays. The gluon system has a smaller 
width in the neutrino channel compared to the hadronic channel where the 
width of the gluon system is broader due to combinatorics. This information 
was not explicitly used in the analysis. Figure [TU] shows the hb, cc and gg 
Higgs mass distribution in the hadronic channel after kinematic fitting. 
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total tagged bs jet 1 after cuts^ 



total tagged cs jet 1 after cuts 




Figure 4: Distributions of flavour tagging variables for bb in hadronic channel: 
Leading jet (a) b-tag; (b) c-tag; (c) c-tag with b background only. 
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I total tagged bs jet 2 after cuts 



I total tagged cs jet 2 after cuts 





Figure 6: Distributions of flavour tagging variables for gg in hadronic channel: 
Leading jet (a) b-tag; (b) c-tag; (c) c-tag with b background only. 
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I total tagged bs jet 2 after cuts 



I total tagged cs jet 2 after cuts 





di-jfit Invariant Mass after classification | 



100 
50 









Mean 115.0± 0.2 




Sigma 7.261±0.161 


LJ 

: , , , , , , . .J, , . \_ 


1 , , , , 1 , , , , 



(a) 



(b) 



di-jet Invariant Mass after classification | 

500 — 




- Mean 115.2± 0.1 



- Sigma 6.450 ±0.096 



(c) 



Figure 8: Di-jet invariant mass (GeV) in the neutrino channel: (a) hh\ (b) cc 
and (c) gg. 
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Mass of higgs2 | 



Mass of higgs2 | 




(a) 

Mass of higgs2 ] 



350 - 
300 3 
250 ; 
200 I 
150^ 

100 ; 

50 = 




Mean 116.4+ 0.2 



Sigma 10.002 + 0.158 



100 150 200 250 



(b) 




Mean 116.2+ 0.1 



Sigma 10.240 + 0.115 



100 150 200 250 



(c) 



Figure 9: Di-jet invariant mass of the Higgs candidate (GeV) in the hadronic 
channel before kinematic fitting: (a) bb; (b) cc and (c) gg. 
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H mass after fitting~| 



H mass after fitting~| 




- Mean 121.32 + 0.05 



- Sigma 8.219 + 0.058 



100 150 



(a) 

H mass after fitting | 



700 E 
600 E 

500 ; 

400 E 
300 E 
200 E 

100 E 



50 100 



200 
100 













Mean 121.5 + 0.2 






Sigma 7.784 + 0.210 


1 


\ 



150 200 



(b) 




- Mean 121.0+0.2 



Sigma 9.961 + 0.237 



200 250 



(c) 



Figure 10: Hadronic channel after kinematic fitting: (a) b di-jet mass; (b) c 
di-jet mass; (c) gluon di-jet mass. 
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2.3 Determination of Branching Ratios 

The branching ratio of the Higgs boson decay to quarks and gluons was 
calculated using events that passed the final neural network selection. The 
calculation was done by normalising the signal cross section to the inclusive 
Higgs cross section, azH = 209±9.8 fb, as determined in an independent recoil 
mass analysis performed for the SiD Letter of Intent [H] . The branching ratio 
is then given by 

BR{H ^ //) = ^ (2) 

(yzH 

where / represents the daughter decay products from the Higgs. The relative 
accuracy of the the branching ratio takes into account both the relative signal 
cross section uncertainty and the relative Higgs-strahlung uncertainty given 
as 

^BR 



BR 

with the relative signal cross section uncertainty calculated by 



AcTj^// signal + background 
(THff signal 

and the cross-section is calculated as follows 

^Hff = Y 

where is the number of signal events after all selections, e is the efficiency 
of signal selection and L is the total integrated luminosity. 

The uncertainty in the efficiency is considered negligible, relying on sim- 
ulations to determine it with sufficient precision. The systematic effects or 
contributions of the luminosity uncertainty were not considered in this anal- 
ysis. 

The weighted average of the signal cross section and its uncertainty are 
calculated using cross section and relative uncertainty values obtained from 
the neutrino and hadronic channels. The weighted average cross section is 
given by 

_ x{6yf + yj^dx)"^ 

Oaverage- ^5^)2 + (5^)2 ' W 

where x and y are the cross sections in the neutrino and hadronics channels 
respectively, and 5x and 5y are the cross section uncertainties in the neutrino 
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and hadronic channels respectively assuming that the two channels are sta- 
tistically independent. The uncertainty of the average cross section is then 
calculated as 

r 6x*dy 

Sz = =, (7) 

where 6z is the uncertainty of the weighted average cross section. 

The events remaining after preselection are categorized using two neural 
networks implemented in FANN [15]. For bb the first NN is trained to dis- 
tinguish the SM background from the inclusive Higgs sample and to produce 
the ^'^sM-Higgs output. In the gg case, the first NN is trained to distinguish 
the SM background from the signal sample and to produce the '^'^sig-SM 
output. The second NN is trained to distinguish the signal from the Higgs 
background sample and to produce the ^^mggs-signai output. The training 
is done separately for bb and for gg, i.e the training is done twice for the first 
NN and twice for the second NN. Figures [TT] and [12] show one-dimensional 
histograms of the trained neutral network outputs. Two-dimensional plots 
of the outputs of the trained NNs are shown in Figures [151 [13 [13 and [T8l 
Again, all the histograms are normalized to luminosity 250 fb^^. 

The final event samples are obtained after applying cuts on First and 
Second NNs. These cuts are optimized by choosing values of the NNs that 
maximise the signal-to-noise ratio as shown in Figures [12] and [TU 

A summary of the results obtained for both the neutrino and hadronic 
channels is given in Tables [6] and [7] 





Neutrino 


Hadronic 


Combined 


Signal events 


2833 


8122 




SM background events 


220 


4700 




Higgs background events 


55 


423 




Signal efficiency % 


24.465±0.004 


26.213±0.002 




Signal amb 


142.7±2.3 fb 


142.5±1.9 fb 


142.57±1.61 fb 


Relative uncertainty on anbb 


1.9% 


1.4% 


1.1% 



Table 6: Measurement results of H— )■ bb branching ratio. 



Table [H] shows the summary of the uncertainties of the Higgs branching 
ratios to bb and gg. Also shown in the table is the uncertainty of the Higgs 
branching ratio to cc as given in |1] 
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Figure 11: bb: Neutrino channel (a) and (b), and Hadronic channel (c) and 
(d). First NN corresponds to (a) and (c), and Second NN corresponds to (b) 
and (d). Solid curves are SM background, dashed curves are inclusive Higgs 
sample (with removed signal) and filled histograms are the signals. 
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I NN signal weFglited~| | NN2 signal weighted | 




(c) (d) 

Figure 12: gg: Neutrino channel (a) and (b), and Hadronic channel (c) and 
(d). First NN corresponds to (a) and (c), and Second NN corresponds to (b) 
and (d). Solid curves are SM background, dashed curves are inclusive Higgs 
sample and filled histograms are the signals. 
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Figure 13: Neural Network Cut optimization for bb : neutrino channel (a) 
and hadronic channel (b) 
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(a) (b) 

Figure 14: Neural Network Cut optimization for gg : neutrino channel (a) 
and liadronic channel (b) 
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"0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



"0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

nn;„.h,„. 



(a) 



(b) 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



(c) 



Figure 15: Neutrino channel bb: Second NN versus first NN for Signal (a) 
Standard Model background (b) and Higgs background (c) 
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(c) 

Figure 16: Hadronic channel bb: Second NN versus first NN for Signal (a) 
Standard Model background (b) and Higgs background (c) 
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rNfi 



0.4 

0.3 |f 




(a) 



NN, 



'Slq-S 





















0.1 0.2 0.3 0.4 0.5 0.6 


0.7 0.8 0.9 



'Siq-SM 



(b) 



|NNhio„.-.,oI 

1p 
















'O 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



(c) 



Figure 17: Neutrino channel gg: Second NN versus first NN for Signal (a), 
Standard Model background (b) and Higgs background (c) 
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Figure 18: Hadronic channel gg: Second NN versus first NN for Signal (a), 
Standard Model background (b) and Higgs background (c) 
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Noutriuo 


Hadronic 


Conil)in(xl 


Signal events 


32 


524 




SM background events 





3621 




Higgs background events 


4 


1431 




Signal efficiency % 


3.245±0.006 


17.673±0.007 




Signal aHgg 


15.1±1.9 fb 


15.6±2.6 fb 


15.41±1.74 fb 


Relative uncertainty on angg 


18.7% 


14.2% 


11.3% 



Table 7: Measurement results of H— T-gg branching ratio. 



Channel 


ABR 0/ 


bb 


4.8 


H— )■ cc 


8.4 


H-^ gg 


12.2 



Table 8: Uncertainties of hadronic Higgs decay branching ratios. 

3 Conclusion 

The measurement of the Higgs boson decay branching ratios to bottom 
quarks and gluons, for a neutral SM Higgs boson of mass 120 GeV, has been 
studied at a centre-of-mass of energy of y/s — 250 GeV and a total integrated 

luminosity of 250 J fb~^. The analysis is based on full detector simulation and 
realistic event reconstruction. The uncertainties on the branching ratios are 
found to be 4.8% and 12.2% for bb and gg respectively. The uncertainty 
in the bb branching ratio is dominated by the uncertainty on the inclusive 
Higgs-strahlung cross section. A good performance of flavour tagging and 
the use of neural networks in event selection are critical in obtaining these 
results. 
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